Requirement of a Primordial Magnetic Field in Chameleon Models 
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We show that the presence of primordial magnetic fields (PMFs) are imperative in order to satisfy 
constraints placed on chameleon models from variation of particle masses at big bang nucleosynthe- 
sis. For initial field values after inflation of order Mpi, the combined magnetic field strength Bo and 
chameleon-photon coupling 1/M F must satisfy (B /5nG) 2 (M F /10 9 GeV) _1 > 2.4 x 1CT 6 . Combin- 
ing these constraints with those derived from considering the degree of mixing between chameleons 
and CMB photons in a PMF, implies the existence of a primordial magnetic field with strength 
Bo > O.lnG, for theories of Modified Gravity with a chameleon mechanism to be viable. 



Gravity theories which extend General Relativity by 
introducing a new degree of freedom have received in- 
creased attention lately due to combined motivation com- 
ing from high-energy physics, cosmology and astrophysics 
PQ . If they claim to account for dark energy or dark mat- 
ter, then they are only valid if they pass the many weak- 
field limit tests. For many this is only possible if they 
have a chameleon mechanism [2]. Such mechanisms 
have the effective mass of the scalar degree of freedom 
being a function of the curvature (or energy density) of 
the local environment, so that in effect the mass is large 
at Solar System and terrestrial curvatures and densities, 
and small at cosmological curvatures and densities. It 
is short-ranged in the Solar System, but becomes long- 
ranged at cosmological densities and can drive the accel- 
eration of the universe. 

In the Einstein frame, the new degree of freedom can 
be identified with the chameleon field. The early universe 
behaviour of the chameleon, considering only the inter- 
action of the chameleon to matter species, was investi- 
gated in [3J. Assuming that the chameleon is produced 
at some phase transition during inflation, and that at the 
end of inflation the chameleon is left at some arbitrary 
field value, it then rolls to the minimum of its potential 
and in some cases oscillates about it. Variation of parti- 
cle masses after big bang nucleosynthesis (BBN) requires 
the maximum amplitude of the chameleon oscillations to 
be below a certain threshold at BBN. Satisfying these 
constraints is of utmost importance in building a viable 
theory of Modified Gravity. In [3J the initial conditions 
of the chameleon field needed to be strongly fine-tuned 
at the end of inflation to satisfy the BBN constraints. 
In this Letter we show that in general the presence of a 
large scale PMF in the early universe is required to drive 
the chameleon towards the minimum in time for BBN. 

The existence of a large scale magnetic field in the early 
universe is a matter of current debate. The WMAP 7- 
year results place a maximum bound on its magnitude of 
Bo < 5.0nG 8 \. Although the presence of a PMF has not 
been detected, it is generally believed to be required for 
the formation of the order fj,G magnetic fields observed in 
galaxies and galaxy clusters. See |13j for a recent review. 



In the chameleon model the matter particles couple 
indirectly to the chameleon <p via a conformal metric 
B^fyg^v. The coupling of the chameleon to charged 
matter naturally generates a direct coupling between 
the chameleon and electromagnetic field of the form 
B F ((j))F^F^ u @]. We take the coupling functions to 
be of generic exponential form: B m (<fi) — exp(</>/M) 
and Bp((j)) ~ exp^/Afp)- The strongest bounds on the 
matter coupling come from particle physics experiments, 
M > 10 4 GeV [5]. Stronger constraints have been derived 
on the chameleon-photon coupling strength. Constraints 
on the production of starlight polarization in the galactic 
magnetic field place a lower bound of Mp > 1.1 x 10 9 GeV 
[6], while measurements of the Sunyaev-Zel'dovich ef- 
fect in galaxy clusters places a lower bound in the range 
M F > (0.25-1.14) x 10 9 GeV, depending on the model 
assumed for the cluster magnetic field [7J. 

The equation describing the chameleon evolution in the 
presence of a magnetic field is, 

- cj>. tt - 3H<P it = V'(4>) + \\BfB' F ^) - T;£C(0), (1) 

where we have assumed a spatially flat FRW metric. The 
universe becomes a good conductor soon after the end of 
inflation, and from then on all magnetic fields are frozen- 
in with |B| 2 oc a -4 . The evolution of the scale factor 
a(t) is described by the Hubble expansion, H = a/a. At 
these early times the universe is radiation dominated and 
we can approximate H 2 (a) ~ H 2 il r oa~' 1 , where fl r o = 
8irGp r o /3Hy is the fractional energy density of radiation. 

In general the stress-energy tensor — — p m , where 
p m is the energy density in non-relativistic matter. How- 
ever as the universe cools and the particle species drop 
out of thermal equilibrium, there is an additional short- 
lived boost to Tff which has a significant effect on the 
chameleon evolution. For the chameleon mechanism we 
require the self-interaction potential of the scalar field 
to be of runaway form. A typical choice is, V(4>) — 
A 4 exp (A™ /4> n ), where n ~ 0(1) and we require A ~ 
2.4 x 10~ 3 eV for a suitable dark energy candidate. 

Equation ([I]) is that of a weakly damped oscillator, 
moving in an effective potential defined by, 

VW(0) = A 4 exp ($) " + !|B| 2 exp (^) +p m exp (£) , 
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where for now we neglect the extra contribution to 
as particle species go non-relativistic. When in equilib- 
rium the chameleon seeks the minimum of this effective 
potential, </> m i n . Notice that M, Mp 3> A and so the 
lefthand side of the potential with <f> < <f> m i n is many 
orders of magnitude steeper than that for <j> > 4>m.in- 
The exponential decay of the functions means that for 
4> < 0min the potential is almost entirely dominated by 
V{4>), while for (j) > 0min the dominant driving force is 
from the p m and \B\ 2 terms. An exact solution to the 
differential equation cannot be found but we use the 
following approximation to determine the evolution of the 
chameleon. Whenever <f> > m ; n we assume V(<j>) is neg- 
ligible. Approximating B m (<f>), Bp(4>) w 1, we can solve 
for the evolution along this side of the potential: 



mum a m ; n is determined by, 



da 

4>{a) 



= -[3a 1 + Aa 2 , 

= -^loga- Aa^ 1 + B, 
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where [3 = \B\ 2 /2M F Q, l0 H$ . We neglect the contribu- 
tion from p m since this is subdominant prior to BBN. 
The constants A and B are determined by the initial 
conditions. Conversely, when <fi < (f> m i n the \B\ 2 term 
is negligible compared to V(<fi). The friction term pro- 
portional to H in equation is also negligible since 
the steepness of the potential dominates. The roll up or 
down on this side of the potential is very rapid and we 
treat it as instantaneous. As the chameleon approaches 
from larger field values, V(<p) effectively acts as a perfect 
clastic collision. 

We assume the chameleon starts from rest at some ini- 
tial value 4>i at the end of inflation. For starting values 
with <pi < 4> m i n , the chameleon falls very quickly to the 
minimum and overshoots to the other side of the poten- 
tial before coming to a halt. It then starts to fall back 
towards the minimum and evolves identically to the case 
with <pi > (j) m in- Depending on the size of the back- 
ground magnetic field and the initial starting value for 
the chameleon field, the chameleon either falls to the 
minimum and start oscillating about it, or the friction 
term dominates and the field remains frozen at its initial 
value. The amplitude of the oscillations gradually decay 
due to the damping term from the Hubble expansion. 

Consider a single oscillation which starts from rest at 
cj> = (f>i at some time a\ with <pi 3> 4>m\n- The field is 
governed by equation ^ as it rolls to the minimum. We 
then assume at 4> m in that it undergoes an instantaneous 
perfect elastic collision and rebounds with the same ve- 
locity. It comes to a halt at the retracement point <p 2 at 
some time a2- 
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where the scale factor when the field reaches the mini- 



log 



Oi 



Oi 



1. 



(5) 



This evolution from one retracement point (</>i,ai) 
to another (0 2 ,02) can be iterated to determine the 
chameleon behaviour. An analytic approximation for the 
evolution of the retracement point (j) max exists in the limit 
of fast oscillations when a m ; n /ai « 1 + S, S <C 1: 
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This power law behaviour is independent of the strength 
of the magnetic field. In this regime, the magnetic field 
strength dictates how rapidly the chameleon oscillates 
but not how long it takes to settle to the minimum. 

In addition to these oscillations there is the extra con- 
tribution to as particle species drop out of thermal 
equilibrium. This is a short-lived boost to the driving 
term on the righthand side of equation 0. We follow 
the convention in [3] and refer to these contributions as 
'kicks'. The first 'kick' occurs at a redshift of approxi- 
mately z ~ 10 14 as the top quark drops out of thermal 
equilibrium. The contribution as each particle, labelled 
by k, goes non-relativistic is [51 111). 
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and the ± sign is for fermions and bosons respectively. 
The t function is approximately 0(1) for x *~ 1 and negli- 
gible otherwise, and so the maximum contribution occurs 
when the temperature of the universe has cooled suffi- 
ciently for it to match the mass of the particle, ~ T. 
The mass of the different species, along with the ratio of 
the number of degrees of freedom to the effective number 
of relativistic degrees of freedom, gk/g*, are in Table m 

An exact solution to the chameleon evolution with the 
'kicks' is not possible. Instead we consider two limiting 
cases: when the duration of the 'kick' is very much less 
that the oscillation period of the chameleon, we can ap- 
proximate the contribution to as a delta function; on 
the other hand when the chameleon is oscillating rapidly 
about the minimum, we can make an adiabatic approxi- 
mation. The duration of the 'kick' can be approximated 
as aija\ rs 100, where a 2 — a\ = Aakkk, for which the r 
function has dropped to less than 10% of its maximum 
value. 

Approximating t(x) as a delta-function centered on 
a = ak, we have for each 'kick', 
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Particle 


mass, rrifc 


gk/g* (at T = m k ) 


type 


t 


173 GeV 


12/106.75 


fermion 


Z 


91 GeV 


3/95.25 


boson 


W ± 


80 GeV 


6/92.25 


boson 


b 


4 GeV 


12/86.25 


fermion 


T 


1.7 GeV 


4/75.75 


fermion 


C 


1.27 GeV 


12/72.25 


fermion 


7T 


0.14 GeV 


3/17.25 


boson 




0.105 GeV 


4/14.25 


fermion 


e 


0.5 MeV 


4/10.75 


fermion 



TABLE I: List of particle species that provide a 'kick' to the 
chameleon evolution as they drop out of thermal equilibrium. 
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X=1, 



1 E+006M and M = 1 E+01 8GeV 



where a k = Tg/rrife, and T 4.9-4.6 is the total area 
under the t(x) curve. We find that the jumps in the 
chameleon field and its velocity from their free evolution 
as a result of the kick are, 



X = 1 E-008, <|>. = 1 E+008M and M = 1 .1 E+009GeV 



A(f> 



-K k M P \Ta k a ' 
-K fc M P1 r (1 - 



where n k = A5M Pl (g k /g+) /tt 4 M. 

When the duration of the 'kick' is significantly greater 
than the oscillation period of the chameleon, we use an 
adiabatic approximation. Treating r(a/a k ) as constant 
over one oscillation, we find that the evolution is identical 
to equation ^ but with /3 replaced by, 



ell 



/3 + M P] J2 K kT 



where a\ is the scale factor at the start of the oscilla- 
tion. In this regime the effect of the kicks is similar to 
a temporary increase in the magnetic field strength, and 
generally boosts the oscillations into the 'fast oscillations' 
regime described by equation 

In Figure [l] we plot two examples of the predicted 
chameleon evolution from the end of inflation (z ~ 10 23 ) 
until BBN (z ~ 10 s ) for different parameter values. Note 
that the evolution depends on /3 and as such is degenerate 
between B$ and M-p. The parameter, 



A= (St) 



Mf V 
1.1 x 10 9 GeV ) 



characterizes this dependence. 

Nucleosynthesis constrains the variation in particle 
masses from BBN until the present day to be less than 
approximately 10% [12]. The particle masses in the 
chameleon model are dependent on <fi with m((f>) — 
B m (cj))mo 0, and so variation in <f> is constrained to be 
less than 0.1M. Since <^ ay) « 10 8 A < 0.1M, this 
leads to 0( BBN ) < QAM. 




FIG. 1: Chameleon evolution for different A, M and </>i. The 
dashed blue line is the evolution of 4>mi n - Green circles mark 
the actual location of the field including the oscillations, while 
the green line is the evolution of </f> ma x in the 'fast oscillations' 
regime. Red crosses mark the location of the 'kicks'. 



We determine the parameter sets that satisfy these 
BBN constraints and plot the exclusion bounds in Fig- 
ure [2] for two choices of matter coupling strength: M = 
1.1 x 10 9 GeV and M = M P \. The ratio of fa to A dictates 
whether the chameleon has started to oscillate prior to 
the top quark dropping out of thermal equilibrium (upper 
plot of Figure [I]) in which case it quickly enters the 'fast 
oscillations' regime with the extra boost from the kicks; 
or whether the field is frozen at its initial value until the 
first kick occurs and sends it towards the minimum. We 
find that in all cases when the field is frozen at its initial 
value, the BBN constraints are not satisfied. The re- 
bound at the minimum of the potential often causes the 
chameleon to shoot to even larger field values by the end 
of the 'kicks' (lower plot in Figure [I]). The field is frozen 
at its initial value for all cases when the interaction with 
the PMF is neglected, and it requires extreme fine tuning 
of the initial conditions to satisfy the BBN constraints. 

The transition between an oscillating chameleon and a 
frozen chameleon prior to the kicks occurs when the field 
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BBN exclusion bounds for M=1 .1x1 GeV and M=M„ 




FIG. 2: The parameter space of 4>i an d A that is excluded by 
BBN constraints. The dashed blue line marks the excluded 
region for M = 1.1 x 10 9 GeV, and the dotted red for M = 
Mpi. The green line marks the transition from the chameleon 
oscillating prior to the 'kicks', and the one in which the field 
is frozen prior to the kicks. 



has oscillated exactly once between a ~ 10~ 23 and the 
first kick at a = 3.8 x 10~ 15 . Under the approximation 
flmin ^ flij we find this transition occurs when, 

(pi « 3.8 x 10 5 Afpi A. 

The transition between when the chameleon is oscillating 
and when it is frozen at its initial value also marks the 
BBN exclusion bound for small but, as we can see 
from Figure [2] there is a saturation point dependent on 
M but independent of A which excludes larger values 
of 4n. This occurs when the chameleon is in the 'fast 
oscillations' regime described by equation ^ for all of 
its evolution and so we have, 



For initial field values up to order Mpi, we find that the 
chameleon to matter coupling strength is constrained at 
M > 10 9 GeV, which is many orders of magnitude greater 
than the existing direct constraints on M. In addition 
there is a degenerate constraint on the PMF strength 
and chameleon-photon coupling, B^qM^ 1 > 2 Ax 1CT 6 . 
This can be compared to the constraints placed on Bq 
and My from considering CMB photons mixing with 
chameleons in a PMF [10]: B 5nG M~ 2 < 6 x 10" 7 at 95% 
confidence (for a PMF spectral index of tib — —2.9). 

The allowed degree of chameleon-photon mixing in the 
early universe places an upper bound on the strength of 
the interaction between the chameleon and a PMF, while 
for the chameleon to be well behaved in the early uni- 
verse and satisfy constraints at BBN we require a strong 
magnetic interaction to drive the chameleon to its mini- 
mum in time. Combining the two bounds, we find there 
is only a narrow band of parameter values that satisfy 
these constraints, 

1.3 x 10 4 B| nG < M 9 < 4 x 10 5 i? 5 2 nG . 

For primordial magnetic fields smaller than 10 _lo G, 
there are no coupling strengths that can satisfy the com- 
bined bound, and so in all cases we require Bq > 10 _10 G. 
This result opens up an altogether new window which 
might lead to a completely different view of the role 
played by early universe magnetic fields when building 
viable extensions to General Relativity. 
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To satisfy the BBN constraints this imposes, <pi < 10 9 M. 
Approximating these bounds as straight lines, we obtain 
the following constraints on M, M-p and Bq for a given 



m > irr 



and B 5 2 nG M 9 - 1 >2.4x 10- 6 -^, 



where S 5nG = S /5nG and M g = M F /10 9 GeV. For 
example when 4n ~ O(Afpi) we require, 

M > 10 9 GeV, and B^JI^ 1 > 2.4 x 1(T 6 . 

In this Letter we have demonstrated that it is essential 
to have a large scale magnetic field in the early universe 
for the chameleon to satisfy bounds placed on variation 
of particle masses at BBN. The results of [3] that did 
not require a PMF are strongly fine tuned to specific ini- 
tial conditions for the chameleon at the end of inflation. 



[1] B. Jain, J. Khoury, Annals Phys. 325, 1479-1516 (2010). 
[2] J. Khoury, A. Weltman, Phys. Rev. Lett. 93, 171104 

(2004); Phys. Rev. D69, 044026 (2004). 
[3] P. Brax et al., Phys. Rev. D70, 123518 (2004). 
[4] P. Brax et al., Phys. Lett. B699, 5-9 (2011). 
[5] D.F. Mota and D.J. Shaw, Phys. Rev. Lett. 97 151102 

(2006); Phys. Rev. D.75, 063501 (2007); P. Brax et al., 

Phys. Rev. D76, 124034 (2007). 
[6] C. Burrage, A. Davis, D. Shaw, Phys. Rev. D79, 044028 

(2009). 

[7] A.C. Davis, C.A.O. Schelpe, D.J. Shaw, Phys. Rev. D80, 

064016 (2009); Phys. Rev. D83, 044006 (2011). 
[8] D. Paoletti, F. Finelli, Phys. Rev. D83, 123533 (2011). 
[9] Y. Fujii and K. Maeda, The Scalar-Tensor Theory of 
Gravitation, Cambridge University Press (2003). 
[10] C.A.O. Schelpe, Phys. Rev. D82, 044033 (2010). 
[11] T. Damour, A. Polyakov, Nucl. Phys. B423, 532 (1994); 

T. Damour, K. Nordtvedt, Phys. Rev. D48, 3436 (1993). 
[12] T. Dent, S. Stern, C. Wetterich, J. Phys. G35, 014005 
(2008). 

[13] R.M. Kulsrud, E.G. Zweibel, Rept. Prog. Phys. 71, 
0046091 (2008). 



